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Abstract. Despite the nice features of the DvaH, Gabadadze and Porrati (DGP) model to 
explain the late-time acceleration of the universe, it suffers from some theoretical problems like 
the ghost issue. We present a way to self-accelerate the normal DGP branch, which is known 
to be free of the ghost problem, by means of an f(R) term on the brane action. We obtain the 
de Sitter self-accelerating solutions of the model and study their stability under homogeneous 
perturbations. 



1. Introduction 

One of the most puzzling problems nowadays in physics is the issue of the late-time acceleration 
of the universe [Ij. A possible approach to tackle this problem is within the frame- work of 
self-accelerating universes [21 [3l IU [5]; i.e. could a modification of gravity at late-time and on 
large scale be the cause of the current inflationary phase of the universe? In other words, could 
this modification of gravity on large scales provide an effective negative pressure that would fuel 
the late-time acceleration of the universe? 

A possible approach to tackle these questions is the Dvali, Gabadadze and Porrati (DGP) 
scenario [2], which corresponds to a five-dimensional (5D) model. In this model, our universe 
is a brane; i.e. a 4D hyper-surface, embedded in a flat space-time. The DGP model has 
two types of solutions: the self-accelerating branch and the normal one. The self-accelerating 
brane is asymptotically de Sitter. This feature takes place without invoking any unknown dark 
energy component. On the other hand, the normal branch requires a dark energy component 
to accommodate the current observations ^j. Despite the nice features of the self-accelerating 
DGP branch, it suffers from serious theoretical problems like the ghost issue In this paper 
we propose a mechanism to self-accelerate the normal branch which is known to be free from 
the ghost issue [6]. This mechanism will be based on a modified Hilbert-Einstein action on the 
brane f5] and the simplest gravitational option is an f{R) term [1]. 

2. The model 

We consider a brane, described by a 4D hyper-surface (/i, metric g), embedded in a 5D bulk 
space-time {B, metric g^^^), whose action is given by 

S= ^d^^yi^ | + ^d^XV=^ j-^ET +a/(ii)+£^| , (1) 

where k| is the 5D gravitational constant, R[g^^^] is the scalar curvature in the bulk and K 
the extrinsic curvature of the brane in the higher dimensional bulk. On the other hand, R is 



the scalar curvature of the induced metric on the brane, and a is a constant that measures 
the strength of the generalised induced gravity term f{R) and has mass square units. Notice 
that therefore the function f{R) has mass square units. Finally, Cm corresponds to the matter 
Lagrangian of the brane. We will assume that the brane splits the bulk in two symmetric 
pieces. The previous action, includes as a particular case the DGP model [2] for f{R) = R and 
a = 1/2k\ where k\ is proportional to the 4D gravitational constant. It can be shown that the 
total energy density of the brane is conserved (we refer the reader to [S] for more details). In 
particular, the energy density of matter on the brane is conserved. 

In what follows, we consider a homogeneous and isotropic brane. The matter sector on the 
brane can be described by a perfect fluid with energy density p^^^ and pressure p^™-), where 
p{m) conserved as we have pointed above. On the other hand, an effective energy density and 
an effective pressure associated to the energy momentum tensor coming from the f{R) term on 
action can be defined as follows [5] 
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p^f^ = 2a|(2ij + 3if2 + ^^/'_i(i?/'_/)[i?/'' + (i?)2;- + 2i/i?r]|, (3) 

Notice that the definition of p^^^ and p^^^ is different from the standard 4D definition in f{R) 
models [5j. On the other hand, the energy density is conserved on the brane. 
The modified Friedmann equation on the brane can be written as 

3^^ = (4) 
While, the spatial component of Einstein equation can be expressed as 

2H + 3H^ + ^ = -^p{p + 2p), (5) 

where the energy density p and the pressure p are defined as 

p = + , p = + p(f) . (6) 
For simplicity, on equations dH) and ([6]) we have used the spatially flat chart of the brane. 



3. Self-accelerating branes 

A de Sitter universe is the simplest cosmological solution that exhibits acceleration and therefore 
it is worthwhile to prove the existence of this solution in our model and study its stability. This 
would be a first step towards describing in a realistic way the late-time acceleration of the 
universe in an f{R) brane- world model. This approach will also enable us to look for self- 
accelerating solutions on the modified normal DGP branch. So, in this section, we first obtain 
the fixed points of the model corresponding to a de Sitter space-time and then we study their 
stability under homogeneous perturbations. 



3.1. Background solutions 

In our model, the Hubble parameter for de Sitter solutions can be expressed a^ 

^ For a maximally symmetric brane in our model, the matter content of the brane behaves like a cosmological 
constant. As such a term can always be reabsorbed in the f{R) term we will disregard the matter content in our 
analysis of de Sitter branes. 



24a^FiHl = 1 + ^Kia'FoiR^Fo - fo) + 1 + '^4aFo [a{RoFo - fo)] (7) 

where e = ±1, the subscript stands for quantities evaluated at the de Sitter space-time, 
Rq = 12Hq and F = df /dR. We recover the DGP model for f{R) = R- In fact, in that case, the 
de Sitter self-accelerating DGP branch is obtained for e = 1 and the normal DGP branch or the 
non-self-accelerating solution for e = — 1. When the brane action contains curvature corrections 
to the Hilbert-Einstein action given by the brane scalar curvature, the branch with e = — 1 is no 
longer flat and accelerates (cf. Fig. [T]). Therefore, an f{R) term on the brane action induce in 
a natural way self-acceleration on the normal branch. Most importantly, it is known that such 
a branch is free from the ghost problem (see [6j and references therein). The reason behind the 
self-acceleration of the generalised normal brane is the presence of the effective energy density 

pi') = a{FoRo - fo) (8) 

on the modified Friedmann equation on the brane. This can be easily shown by comparing the 
Friedmann equation ([7]) with that of modified gravity on brane world-models [7] 
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Figure 1. The figure on the left shows the behaviour of the rescaled squared Hubble rate 
2Kfa'^FQHQ for the two branches that generalise the DGP solution versus the rescaled energy 
density p^'^' defined Fq{RqFo — fo). The blue star corresponds to the normal DGP 

branch which is flat. The red star corresponds to the self-accelerating DGP branch. On the 
other hand, the blue curve corresponds to the generalised (by the inclusion of the f{R) term) 
self-accelerating branch, while the red curve corresponds to the generalised (by the inclusion of 
the f{R) term) normal branch. The figure on the right corresponds to a zoom of the normal 
branch as it appears on the figure of the left. 



3.2. Stability of the self- accelerating solutions 

We next analyse the stability of de Sitter solutions under homogeneous perturbations up to first 
order on 5H = H(t) — Ho- We will follow the method used in [8]. 

The perturbed Friedmann equation @ implies an evolution equation for 5H: 



6H + 3Ho6H + mlg6H = 0, 



(9) 



where m^g is defined as 



where 
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and /rr = f /dB?. All this quantities are evaluated at the de Sitter background solution. Any 
de Sitter solution is stable as long as is positive. 



The terms defined on Eq. (jlip have the following physical meaning: (i) 
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is the analogous 



quantity to m^g in a 4D f(R) model [8] , (ii) m-back ^ purely background effect due to the shift 
on the Hubble parameter respect to the standard 4D Case and (iii) rn^^^^ is a purely perturbative 
extra-dimensional effect. 

If we assume that we are close to the 4D regime; i.e. the Hubble rate of the brane is close to 
its analogous quantity in a 4D f{R) model, then rn^^y^ > and rrip^^^ < 0. Consequently, rn^^^^ 
tends to make the perturbation heavier. However, the perturbative effect encoded on rrip^^^ 
would make the perturbation lighter. It can be shown that the extra-dimension has a benigner 



effect in the 4D f(R) model; i.e. 
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4. Conclusions 

We have presented a mechanism to self-accelerate the normal DGP branch which unlike the 
original self-accelerating DGP branch is known to be free of the ghost problem. The mechanism 
is based in including curvature modifications on the brane action. For simplicity, we choose 
those terms to correspond to an f{R) contribution, which in addition is known to be the only 
higher order gravity theories that avoid the so called Ostrogradski instability in 4D models. 
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We have assumed the natural condition Fq > 0; i.e. the effective gravitational constant of the brane is positive. 
On the other hand, we have also assumed that we are slightly perturbing the Hilbert-Einstein action of the brane, 
i.e. /o ~ Rq. Therefore, fo is positive because i?o 
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